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Abstract. We describe Very Large Array observations of the extreme carbon star IRC-f 10216 at 8.4, 14.9, and 
22.5 GHz made over a two year period. We find possible variability correlated with the infrared phase and a cm- 
to sub-millimeter wavelength spectral index very close to 2. The variability, observed flux densities, and upper 
limit on the size are consistent with the emission arising from the stellar photosphere or a slightly larger radio 
photosphere. 
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1. Introduction 

The "extreme" carbon star IRC+10216 (CW Leonis) was 
disco vered in the early d ays of modern infrared astron- 
omy ijBecklin et al.lll969ll . The star is one of the bright- 
est near- and mid-infrared sources (LeBertre 1987 and 
references therein). IRC+10216 has an extremely rich 
molecular spectrum arising fron i a dense envelope cre - 
ated by a powerful mass outflow ijCernicharo et al1l200Cl() . 
These attributes, together with its proximity (distance 
of « 130 pc) and its prodigious mass-loss rate (1.5 - 
3.3 X 10-^ Mq yr-i), make IRC+10216 a "keystone" ob- 
ject for which many aspects of (carbon-star) asymptotic 
giant branch (AGB) evolution can be studied in detail. 
IRC+10216's properties are summarized in Tabled 

IRC+1 0216 has b een detected at wavelength s of 1.5 
and 2 cm l|Sahai et al. (1989). Drake et all l)l99l[l '). From 
Very Large Array (VLA) observations it should be possi- 
ble to measure its size, position (which is extremely dif- 
ficult at other wavelengths owing to absorption and scat- 
tering) , proper motion, and the spectral index of its radio 
emission. Very Large Array data can strongly constrain 
our information on the star's nature. 

This paper has the following structure: In Sj^J we 
describe six-epoch, three- wavelength VLA observations. 
These, and an additional dataset (published by Drake et 



al. 1991), were retrieved from the VLA archive. The prop- 
erties of the continuum emission are discussed in In 
particular, we address the variability, spectral index, ab- 
solute position, and source morphology. In 21^6 synthe- 
size these results to reach conclusions on the nature of the 
radio emission. 



2. VLA Observations and data reduction 

2.1. 1991 - 1993 observations 

The data reported here were taken by two of us (RS & 
MC) with the NRAO Very Large Array (VLA)^ We an- 
alyzed observations made between 1991 June 1 and 1993 
June 9, spaced typically by 3 to 4 months. Observations 
were made at 8.4, 14.9, and 22.5 GHz (X-, U-, and K- 
band, respectively)^. Details are given in Table 12 A typ- 
ical observing run consisted two ~ 20 - 28 minute long 
X- and U-band observations, preceded and followed by 
observations of the calibrator 0953+254. Because of the 
generally short coherence times at the highest frequency, 
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Table 1. IRC+10216 - Fundamental Parameters 



Name 




R6f6rGIlC6 


Aliases 


CW Leonis, RAFGL 1381 


SIMBAD" 




IRAS 09452+1330 




QJ2000 


09''47™57.'392' 




'5j2000 


+13°16'43"63'' 




Spectral type 


C 




Effective temperature, T^g 


2000 K'^ 


Groenewegen (1997) 


Luminosity, 


8.2 10^ [D{kpc)f Lc/ 


Groeneweeen ('1997") 


Mass loss rate, M 


1.5 - 3.3 X 10"^ A/0 yr~^'= 


Crosas & Menten (1997"). 






^roenewegen et al. (199^ 


Distance, D"^ 


110 - 150 pc 


I2rosas & Menten (1991), 




16 km s"^ 


Groeneweeen et al. (1998) 


Expansion velocity 


Loud et al. (1993) 


LSR velocity, «lsr 


-22 + 1.5 km s"^ 


Loud et al. (1993) 


Heliocentric velocity, «hci 


-15 km s-i 





° |http: / / sirobad^u^strasbg.fr/l^ osition derived from 1993 January 26 VLA data. The estimated uncertainty is 25 milli arcseconds 
(see tliM . ^Groenewegenl (119971) calculates a range of models resulting from a range of values for this parameter; we adopt 2000 K, 
which is used in most of his models. Assuming TcS — 2000 K. "^The listed mass loss rate and distance estimates represent the 
range of values derived from modeling of multi-line CO observations described in the referenced studies. 



the K-band scans were shorter (« 15 minute duration) 
and more scans were made in this band. At the end of 
each run, we observed the absolute flux density calibrator 
3C 286. 

The flux density of 0953+254 was determined by com- 
paring its visibility amplitude with that of 3C 286. The 
elevations of 3C 286 were comparable to that of 0953+254, 
minimizing systematic amplitude calibration errors due to 
elevation-dependent antenna gain. By comparing the vari- 
ations of the flux densities for several scans of 0953+254, 
we estimate our flux density scale is accurate to within 
10% at X-band, 15% at U-band, and 20% at K-band. Data 
from one epoch (1992 October 2) suffered serious instru- 
mental difficulties and were discarded. Table [5] contains 
the measured flux densities of 0953+254. Complex gain 
correction factors were interpolated from the amplitude 
and phase of 0953+254 to the times of the IRC+10216 
data. 

We made maps using the AIPS task IMAGR and 
fit elliptical Gaussian models to the images. The in- 
tegrated flux densities are listed in Table HJ Because 
IRC+10216 is weak, self calibration was not possible, and 
the "fast-switching" technique, routinely employed now, 
had not been developed at the time of our observations. 
Nevertheless, during one epoch (1993 January 26) when 
the array was in the A-conflguration, superb phase coher- 
ence persisted throughout the observations. We use those 
data to establish the source position (in H3.3(l . 

2.2. Other VLA data 

VLA A-configur ation observations of IRC+10216 have 
been reported bv lPrake et alJ ( 199lh . We re- analyzed their 
1987 June 2 and 3 U-band data to compare with our U- 



band data. Fitting an elliptical Gaussian using the AIPS 
task JMFIT, we find the source to be nearly unresolved, 
with upper limits on the deconvolved size of 95 mil- 
liarcesconds (mas) and a peak brightness and integrated 
flux density of 1.29+0.15 mJy beam"! and 1.40+0.27 mJy, 
respectively. The latter value can be compared with our 
observed range of integrated U-band flux densities of 1.77 
to 2.67 mJy (see TableEJ. Since we used identical calibra- 
tion procedures for both datasets and the same absolute 
flux calibrator, this discrepancy is probably real. Since the 
1987 data employed a secondary gain calibrator nearer to 
IRC+10216 than used for our data, one would expect that 
the 1987 data should yield a greater flux density than the 
1993 data, if random atmospheric phase fluctuations de- 
grade the imaging. 

In order to estimate of the absolute position accuracy, 
we made images from the the June 2 and 3 data separately. 
The positions of the fltted Gaussians differed by 18 mas 
in right ascension and 10 mas in declination, testifying to 
the excellent phase coherence during these observations. 
Assuming that large-scale atmospheric effects, which usu- 
ally limit astrometric accuracy, were independent for these 
two days, this indicates an uncertainty of « 20 mas in each 
coordinate. 



3. Properties of the radio emission 

3.1. Variability 

In Fig. n] we show the flux densities of IRC+10216 listed 
in Table 2. The error bars include contributions from ran- 
dom noise and systematic calibration uncertainties. The 
data show probable variability. We flt the flux densities 
of IRC+10216 with a model that includes a constant and 
a sinusoidally varying component. Rather than fitting the 
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Table 2. Log of IRC+10216 multiband VLA observations 



JD 2440000+ 


Civil Date 


Array 


5(60953+254) (Jy) 




5'(IRC+10216)(mJy) 










X-band 


U-band 


K-band 


X-band 


U-band 


K-band 


8410 


1991 June 01 


D(-A) 


1.80(0.01) 


1.90(0.05) 


2.27(0.07) 


0.92(0.07) 


2.67(0.39) 


7.10(0.68) 


8524 


1991 September 24 


A^B 


1.80(0.08) 


2.02(0.13) 


2.40(0.13) 


0.48(0.07) 


2.11(0.55) 


3.87(1.20) 


8652 


1992 February 09 


CnB 


1.56(0.01) 


1.58(0.01) 


1.51(0.01) 


0.82(0.04) 


2.23(0.17) 


4.24(0.33) 


8772 


1992 May 30 


C-^D 


1.54(0.01) 


1.52(0.02) 


1.33(0.03) 


0.86(0.07) 


2.42(0.20) 


3.18(0.67) 


9014 


1993 January 26 


A 


1.58(0.03) 


1.55(0.04) 


1.60(0.05) 


0.60(0.04) 


2.00(0.31) 


5.33(1.03) 


9143 


1993 June 09 


CnB 


1.72(0.01) 


1.69(0.05) 


1.77(0.24) 


0.48(0.04) 


1.77(0.21) 


3.86(0.60) 


Average 












0.67(0.03) 


2.14(0.18) 


4.20(0.46) 



The first six lines present the resuhs of our VLA measurements epoch by epoch. The first and second columns give the JuUan 
and civil dates, respectively. The third column gives the array configuration. Two letters connected by an arrow denote epochs 
when the array was in transition from one to another configuration; CnB denotes a mixed C- and D-configuration. The forth 
to sixth columns give the bootstrapped flux density of the phase calibrator B0953+25 at X-, U-, and K-band respectively, 
while the seventh to ninth columns give the same for IRC+10216. The quoted uncertainties are formal errors from JMFIT. We 
estimate additional systematic errors of 10, 15, and 20% for X-, U-, and K-band data, respectively. Based on the formulae given 
bv Baars et al.l (119771) . the absolute flux calibrator 3C 286 was assumed to have flux densities of 5.21, 3.46, and 2.52 Jy at X-, 
U-, and K-band. For each band, data from two 50 MHz wide channels were taken in right- and left-hand circular polarization. 
The center frequencies were 8.439, 14.940, and 22.460 GHz for the X-, U-, and K-band data. For the first epoch (1991 June 1, 
JD 2448410) only visibilities within a u, ?;-range of 50 kA, corresponding to D-array antenna separations, were considered. Data 
from one more epoch (1992 October 2, JD 2448897) had to be discarded due to instrumental problems. The bottom line gives 
weighted averages, with individual-epoch uncertainties estimated by adding quadratically the formal uncertainties from JMFIT 
with estimated absolute calibration uncertainties. 



data from each band separately, we assume that a constant 
spectral index, a, describes the difference in flux density 
with observing frequency. Specifically, the model for the 
flux density, S{h', t), at a frequency, v, and time, t, is given 
by the following: 

S{iy, t) = (5(i.o) + A5(i.o) cos (At/P)) [v/v^Y , 

where S^vq) and AS^vq) are the constant and amplitude 
of the time varying flux density at a reference frequency 
vq (set to 10 GHz), At is the time offset from maximum 
(imax), and P is the period. 

A weighted least-squares fit to the data yielded the fol- 
lowing results: S'(z^o) = 1.02 + 0.14 mJy, AS'(i'o) = 0.25 + 
0.09 mJy, P = 535 + 51 d, W = (JD2440000+)8841 + 
24 d, and a = 1.88 + 0.14. This fit has a reduced per 
degree of freedom of 1.6, and the parameter uncertainties 
are the formal la values scaled by -\/L6. Thus, IRC+10216 
appears to display flux density variations of +25% at a sig- 
nificance level of nearly 3a. The variability will be further 
discussed in the Appendix. 

IRC+10216 sho ws strong periodic variability at in - 
frared wavelengths l|Dvck et all Il99lt iLe Bertrel IT99I . 
Dyck et al. determine a period of 638 d at (infrared) K- 
band, while Le Bertre's measurements in the J-, H-, K-, 
L'-, and Af -bands yield values between 636 and 670 d. If 
we assume the infrared period has a 5% uncertainty, then 
the cm-wave period of 535 + 51 d differs from the infrared 
period by about 2a. The cm- wave maximum occurs at 
8841 + 24, whereas the infrared maximum (extrapolated 
by about 1200 d) would be predicted to occur near 8770 
with an estimated uncertainty of about +50. Thus, the 
cm-wave minus infrared maxima differ by 71 + 71 d and 
are consistent within their joint uncertainty. 



3.2. Spectral index 

As discussed above, the radio emission from IRC+10216 
is probably variable. Therefore, we calculate the X-U- 
K-band spectral index (SI) for each epoch and flnd val- 
ues between 1.6 and 2.2, with uncertainties of up to 20%. 
These daily Sis are consistent with a constant value, and 
a variance-weighted average of these indices is 1.92 + 0.11. 
We also listed in TableElvariance- weighted, time-averaged 
flux densities for each band. Using these to determine a 
SI yields 1.91 + 0.11. The close correspondence between 
the daily and time averaged Sis suggests that the spectral 
index is not strongly variable. 

The cm-wave spectral index obtained by least-squares 
fltting of all of the data simultaneously with the model 
described above is 1.88 + 0.14. This value is consistent 
with those observed for Mira va riables and with t he valu e 
of 1.86 for a model of Miras iReid fc MentenI l|l997a|) . 
However, the cm-wave spectral index is also consistent 
with a value of 2.0 for optically-thick black-body emis- 
sion. 

At millimeter- and sub-millimeter wavelengths, flux 
densities measured with single-dish telescopes (with > 10" 
FWHM beams) contain contributions from the star as well 
as from the extended dust envelope. If the measurements 
are made with wideband bolometer detectors, circumstel- 
lar molecular emissi on may also contribute » 30% to 
the total flux density ijWalmslev et al.ll99lllGroenewegenl 
ilS97) . Whereas emission by dust dominates in the submil- 
limeter region on large scales (see Fig. 4), it is negligible 
at cm- wavelengths. Based on the model discussed in the 
Appendix, we expect that the dust envelope contributes 



4 



Menten et al.: Radio continuum observations of IRC+10216 



only ^0.1 mjy to the total flux density of ^2 niJy in our 
VLA K-band measurements. 

A few interferometric observations exist that should 
deliver accurate measurements of compact emission free 
of dust contamination. Also by selecting portions of the 
band largely free of molecul ar lines, nearly pure con tinuum 
emission can be measured. iLucas &: GuelinI ()l99 9^ , using 
the IRAM Plateau de Bure Interferometer (PdBI), find a 
"point source" at 89 and 242 GHz with flux densities of 
65 ± 7 and 487 ± 70 mJy, respectively. (Note these uncer- 
tainties are much larger than the formal errors given in 
the cited paper; R. Lucas personal communication.) These 
authors also detect "extended" emission with flux densi- 
ties of 12 and 57% of the point source values at 89 and 
242 GHz, respectively, whic h they ascribe to d ust emis- 
sion from the inner envelope. lYoung et al.l ll2 004l) u sed the 
Smithsonian Sub-Millimeter Arrav ijHo et al.ll200^ to im- 
age the 680 GHz continuum emission of IRC-l-10216. They 
also find a "compact, unresolved component" in their « 2" 
FWHM beam with a flux density of 3.9 ±1.2 Jy. The flux 
densities discussed above are plotted in Fig. |21 A least- 
squares fit to all the data (cm through sub-millimeter) 
delivers an SI of 1.96 ± 0.04. 

3.3. Absolute position of the radio emission 
distribution 

During our 1993 January 26 observations, the VLA was in 
its most extended configuration and the phase coherence 
were excellent. Data from that epoch can be used to ob- 
tain meaningful images of the emission distribution and to 
make an astrometric position determination. To check the 
validity of the size information and the position accuracy, 
we determined, for the K-band data, phase and amplitude 
calibration of only the first scan of 0953-1-254 and then ap- 
plied extrapolations of these solutions to subsequent scans 
of that source, which were spaced by « 16 minutes. As ex- 
pected, the phase coherence degraded with elapsed time 
between the first and subsequent scans. However, in maps 
made from the uw-data of the second scan, taken 16 min- 
utes after the first, the source showed < 3% amplitude 
degradation, a (formal) source size < 31 mas and an offset 
of < 5 mas in both coordinates from the nominal position. 
Since IRC-l-10216 scans were placed between 0953-1-254 
scans, we conclude that some, but probably little, blur- 
ring occurs due to deficient phase calibration. The test 
described above does not account for phase calibration er- 
rors caused by the fact that 0953-1-254 is > 12° away from 
IRG-l-10216. To take that into account, we increase our 
value for the systematic position uncertainty to 20 mas. 

The position uncertainty may be dominated by statis- 
tical errors due to the modest signal to noise ratio of our 
continuum images. JMFIT delivered positions for the X-, 
U-, and K-band images differing by maximally 20 mas 
from each other, with formal errors of between 5 and 
16 mas. The total uncertainty, obtained by quadratically 
adding the systematic and statistical uncertainty is 15 mas 



in each coordinate. The variance-weighted mean position 
and its uncertainty is g iven in Table [D Th e position de- 
rived from the data of IPrake et al.l l|l99lh is (a, <5)j2ooo 
= 09^47™57.''382, +13°16'43"61, for which we estimate an 
uncertainty of 20 mas in each coordinate (see il2.2|l . The 
latter position (measured on 1987 June 2.5) is offset from 
the 1993 January 26 position by (—146,-20) mas in the 
(East, North) directions. Taking these offsets at face value, 
we derive proper motion components of IRC-l-10216 of 
(-1-26 ±6,-1-4 ±6) mas yr~^, respectively. A motion of that 
magnitude is plausible. For example, Mira (o Ceti), which 
is at a similar distance as IRC+10216, has a larger proper 
motion of (+33, - 239) mas yr^^ measured by Hipparcos 
(|Perrvman et alJ lT997) . 

The proper motion derived abo ve is consistent with 
the upper limit of 30 mas yr" -1 that iBecklin et all l|l969tl 
derive from a comparison of a Palomar Sky Survey plate 
taken in 1954 with a plate they took in 1969. 

3.4. IVIorphology 

Our best quality radio data (1993 Jan 26) suggest an elon- 
gated flux distribution (position angle, P. A., 82±5 degrees; 
i.e. roughly east- west direction), with major and minor 
axes of 213 ± 20 and 73 ± 20 mas, respectively, which is 
significantly larger than the 1987 upper limits of 95 mas 
discussed in 42. 2\ 

While an elongated source would be tantalizing, we 
have greater confidence in the 1987 than the 1993 data. 
The separation between IRC+10216 and 0952+179, the 
phase calibrator used by Drake et al., is only 1?7 in right 
ascension and 4?4 in declination. In contrast, the phase 
calibrator used in 1993 is displaced by 2?1 in right as- 
cension and a much larger distance, — 12?2, in declination 
from IRC+10216. While it is difficult to tell conclusively, 
we believe that the source size observed in 1993 January 
may be affected by uncorrected residual phase errors stem- 
ming from the large angular separation of IRC+10216 and 
the calibrator. 

ICroenewegenl lll996l) takes the elongation found by 
IPrake et al.l l|l99lf) at face value and argues for as- 
phericity. We note that on scales of a few stellar 
radii, IRC+10216's morphology is highly asymmetric and 
changes with time. Multi-epoch 2.2 fun speckle interferom- 
etry between 1995 and 1998 shows an asymmetric clumpy 
shell with two dominant components, A and B, whose sep- 
aration grows from 19 1 to 265 m illiarcscconds (mas) be- 
tween 1995 and 1998 (,Osterbart et al. 2000) . 2.2 ^m Keck 
I aperture masking interferometry be tween 1997 and 199 9 
confirm the observed morphology (Tvithill et aLl l200fl|) . 
While both studies find motions in the material, these are 
inco nsistent with simple sph erical expansion and a bipolar 
flow. l Osterbart et alJ l)2000l) suggest that the star itself is 
identic al to, or located nea r, component B. On the other 
hand, l)Richichi et alJl2003|) . modeling their infrared lunar 
occultation data, favor an identification of the star with 
component A. 
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Future VLA observations outlined in ^ together with 
absolute infrared astrometry, may unambiguously estab- 
lish the star's position relative to the clumpy components 
of its envelope. Also, future VLA array observations at 
1.3 and 0.7 cm can definitively settle the question of 
IRC+10216's radio source size and morphology. 



4. The nature of the radio emission 



iGroeneweeenl l)l997l) performed model calculations of 
IRC-|-10216's infrared through radio continuum emission. 
His models and assumptions under-predict our measured 
radio flux densities by a factor of ~ 1.8. A nearer dis- 
tance of « 100 pc (instead of the 135 pc assumed) would 
provide a better fit to the data. For X- and U-band, his 
models predict the emission to arise solely from the stellar 
photosphere, and at K-band the dust contribution from 
the whole envelope would add about 25% to the photo- 
spheric emission. Assuming an effective temperature of 
2000 K, a photosphcric diameter of 70.2 mas is derived, 
which is consistent with our upper limit. As shown in the 
Appendix, assuming these values for temperature and di- 
ameter, IRC-|-10216's spectral energy distribution can be 
modeled successfully. 

A basic question regarding the cm-wave emission is 
whether it comes from the stellar photosphere or a some- 
what larger radio photosphere. The radio emission from 
Mira variables comes from a radio photosphere, which is 
controlled by H~ free-free opacity and is approximately 
twice the size of the stellar phot osphere (defined in line - 
free regions of the IR spectrum) l|R,eid fc Mentenlll997ajl . 
If a size (or stricter upper limit) could be measured with 
the VLA, that would directly yield a brightness tempera- 
ture. A temperature below 2000 K would be expected for 
a radio photosphere, whereas a temperature in excess of 
2000 would be expected for the stellar photosphere. 

The amplitude of the radio variations of IRC+10216 
of ±25% is larger than for Mir a variables, whe r e varia - 
tions of < ±15% are observed bv lReid fc MentenI l|l997aj) . 
This is in contrast to variations of ±40% (see appendix) 
that would be expected from the stellar photosphere of 
IRC±10216, based on IR variations of nearly 2 magni- 
tudes. For Miras, simple models of radio photospheres 
with shocks propagating outward with speeds less than 
7 km s^^ can explain the absence of strong radio vari- 
ability. However, shocks propagating outward with speeds 
of ~ 10 — 15 km s~^ can pr oduce variation of the m agni- 
tude we find in IRC+10216 l|R,eid & Mentenlll997b^ . 



5. Conclusions and outlook 

Using the VLA, observations of IRC+10216 at 8.4, 14.9, 
and 22.5 GHz were made over a three year period. We 
find probable variability and a spectral index of about 
1.9, i.e. near that of optically thick blackbody emission, 
but consistent with a radio photosphere as observed in 
Mira variables. 



Future VLA A-configurations observations at Q-band 
(7 mm wavelength) with a resolution of 35 mas should def- 
initely resolve the emission and yield a measurement of the 
brightness temperature. This could discriminate between 
a stellar and radio photosphere. Such observations may 
even lead to the detection of surfa ce structure, as observed 
for a Orionis bv lLim et al.l l)l998() . These observations will 
benefit greatly by employing the "fast-switching" tech- 
nique. 

Appendix - Modeling IRC+10216's Spectral 
Energy Distribution 

The radiative transfer in the dusty circumstellar envelope 
around IRC+ 10216 has been c o mputed by var i ous au thors 
(for instance. lLe Bertr'3 1 19871) : lGroenewegenl 1 19971) ). but 
th e most ambitious and compl ete modeling was performed 
by iMen'shchikov et al.l l)200l[l who a ) considered a two- 
dimensional configuration to explain the non-radial struc- 
ture close to the star (d iscovered at 2.2 um w avelength 
speckle interferometry bv lOsterbart et all 1 20001) ). b) took 
into account details in the composition of the dust by 
adding, to the major component of amorphous carbon, ad- 
mixtures of silicon carbide and magnesium sulfide grains 
to reproduce broad spectral features seen by the Infrared 
Space Satellite (ISO) and from the ground (at 11 and 
'^SOfim) and, cj, included the phase in the stellar cycle. 
The bolometric luminosity of the star oscillates in their 
model by a factor of 2.5, although the average (infrared) 
if -band variations of '-^2mag (Dyck et al. 1991) would 
suggest a greater ratio of maximum to minimum luminos- 
ity) imax/ ^min — 6. 

We add to the existing radiative transfer calculations 
another one, on a much more modest level, which has the 
goal of clarifying the relative contributions of the star and 
its dusty envelope to the long wavelength emission. In light 
of the huge amount of observational data, simplicity has, 
despite its obvious shortcomings, the advantage of trans- 
parency and of not being entangled by occasionally sec- 
ondary details. The numerical code is based on a standard 
ray tracing method assuming radial symmetry. It can han- 
dle an arbitrary number of grain types, includes scattering 
and is described in section 13.2 of .Kriigel 1,200^ . 

Our model parameters are as follows: The star has 
at maximum a luminosity L = 1.1 x 10^ Lq and emits 
as a blackbody of temperature = 2000 K. The den- 
sity distribution of the dust is radially symmetric with 
p{r) (X . The outer radius of the circumstellar envelope 
equals 3 x 10^^ cm, corresponding to a source diameter of 
6' at a distance of 110 pc, the inner radius is 3 x 10^** 
cm, or about five stellar radii. There the dust tempera- 
ture is around 1000 K, which is close to the expected con- 
densation temperature. The optical depth in the T^-band 
toward the star amounts to Ay — 48 mag. Neglecting mi- 
nor mineralogical components, the dust consists only of 
spherical amorphous carbon grains with a n{a) c 



-3.5 



size distribution and lower and upper limits a_ = 
a+ = 640A, respectively, and optical constants as com- 
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piled by IZubko et all l)l996|) for their type BE. At 1mm 
wavelength, the mass absorption coefficient, A'^, is about 
0.5 cm^ per g of dust with a v^'^ frequency dependence. 
As the grains are much smaller than the wavelength, their 
size distribution does not affect the value of K^, at FIR or 
longer wavelengths. Being a radiative transfer model for 
a dusty medium, the dust-to-gas ratio does not enter our 
calculations. 

Figure 121 displays model results together with observa- 
tional data appropriate near maximum luminosity. The 
solid curve is our fit near maximum luminosity within 
a 6' beam. It was obtained by adjusting the total opti- 
cal depth, the inner radius, and the grain sizes, the other 
above mentioned parameters were kept fixed. Overall, the 
fit is satisfactory. The squares which lie significantly below 
the model refer to smaller angular sizes (see the caption 
of Fig.EJ). 

Figure^ shows how the flux increases with beam size. 
As the curves are fiat in the very left part of the figure, 
the emission from diameters less than 0."2 is due to the 
star, not to dust, despite the central dust density peak. 
One can read from the figure that the fraction of the total 
flux contributed by the stellar photosphere increases from 
-4% at lOO^fm to -96% at 1.3 cm. 

One can qualitatively understand why the bright- 
ness variations over a pulsational period are considerably 
smaller in the radio regime (Fig^l than at near infrared 
wavelengths. The monochromatic flux density which we 
would receive without foreground extinction, i^^, from a 
star with a blackbody photosphere of temperature T and 
luminosity L is 

LB,{T) 



where a denotes the Stefan-Boltzmann constant and D 
the distance. The equation becomes obvious by integrat- 
ing over frequency which yields AttD^F = L, with F = 
J F„dv. In the Rayleigh- Jeans limit, at radio wavelengths, 
B^{T) cx T and therefore F^ L/T^. 

How do L and T vary during a cycle? Reid & Menten 
(1997a) in their Eg . (1), w hich is based on observations by 
IPettit fc Nicholson! ( 1235^ , approximate the stellar radius 
and temperature of a Mira variable as a function of phase 
(/) by T{(j)) = 2300 + 300 cos K and i?(0) = 2 + 0.4 sin 
AU. Oscillations in luminosity, L{<j)), follow because L — 
AnR^aT*. According to these formulae, temperature and 
luminosity are roughly in phase, and therefore the ob- 
served radio flux Fi, (x L/T^ changes less than L itself. 
The radio F^ varies from maximum to minimum by a fac- 
tor of 2 (corresponding to excursions from the mean by 
±40%) and L varies by a factor 3.76. We note that at 
infrared wavelengths, assuming a constant optical depth 
over the cycle and no dust emission, the monochromatic 
flux mimics L almost closely. The dotted curve in Fig|31 
is our model for IRC+10216 at minimum (2.5 times less 
luminous, T = 1700 K). At minimum, the 1.3 cm radio 
brightness is 25% lower and the J- and i7-band flux den- 
sities are 2.5 times lower compared to maximum light. 



The radio variations around the mean suggested by 
Fig^are smaller than ±40%. This suggests that the cm- 
wave emission may not come directly from the stellar pho- 
tosphere. However, other possible reasons for this discrep- 
ancy might be a smaller variability in the bolometric lu- 
minosity (by a factor 2.5, as Men'shchikov et al. (2001) 
assume, and not by 3.76) or inaccuracies in the simple 
expressions for R{4>) and T(0). 

Acknowledgements. We would like to thank the referee for 
comments and, in particular, Malcolm Walmsley for taking the 
editor's job very seriously, which led to a significant improve- 
ment of the paper. This research has made use of the SIMBAD 
database, operated at CDS, Strasbourg, France. 

References 

Baars, J. W. M., Genzel, R., Pauliny-Toth, I. I. K., & Witzel, 

A. 1977, A&A, 61, 99 
Becklin, E. E., Frogel, J. A., Hyland, A. R., Kristian, J., & 

Neugebauer, G. 1969, ApJ, 158, L133 
Cernicharo, J., Guelin, M., & Kahane, C. 2000, A&AS, 142, 

181 

Crosas, M. & Menten, K. M. 1997, ApJ, 483, 913 
Drake, S. A., Linsky, J. L., Judge, P. G., & Efitzur, M. 1991, 
AJ, 101, 230 

Dyck, H. M., Benson, J. A., Howell, R. R., Joyce, R. R., & 

Leinert, C. 1991, AJ, 102, 200 
Groenewegen, M. A. T. 1996, A&A, 305, L61 
— . 1997, A&A, 317, 503 

Groenewegen, M. A. T., van der Veen, W. E. C. J., & 

Matthews, H. E. 1998, A&A, 338, 491 
Ho, P. T. P., Moran, J. M., & Lo, K. Y. 2004, ApJ, 616, LI 
Kriigel, E. 2003, The physics of interstellar dust (loP Series in 

astronomy and astrophysics, ISBN 0750308613. Bristol, UK: 

The Institute of Physics, 2003.) 
Le Bertre, T. 1987, A&A, 176, 107 
— . 1992, A&AS, 94, 377 

Lim, J., Carilh, C. L., White, S. M., Beasley, A. J., & Marson, 

R. G. 1998, Nature, 392, 575 
Loup, C, Forveille, T., Omont, A., & Paul, J. F. 1993, A&AS, 

99, 291 

Lucas, R. & Guelin, M. 1999, in lAU Symp. 191: Asymptotic 

Giant Branch Stars, 305 
Men'shchikov, A. B., Balega, Y., Blocker, T., Osterbart, R., & 

Weigelt, G. 2001, A&A, 368, 497 
Osterbart, R., Balega, Y. Y., Blocker, T., Men'shchikov, A. B., 

& Weigelt, G. 2000, A&A, 357, 169 
Ferryman, M. A. C., Lindegren, L., Kovalevsky, J., et al. 1997, 

A&A, 323, L49 
Pettit, E. & Nicholson, S. B. 1933, ApJ, 78, 320 
Reid, M. J. & Menten, K. M. 1997a, ApJ, 476, 327 
— . 1997b, Ap&SS, 251, 41 

Richichi, A., Chandrasekhar, T., & Leinert, C. 2003, New 

Astronomy, 8, 507 
Sahai, R., Claussen, M. J., & Masson, C. R. 1989, A&A, 220, 

92 

TuthiU, P. G., Monnier, J. D., Danchi, W. C, & Lopez, B. 

2000, ApJ, 543, 284 
Walmsley, C. M., Chini, R., Kreysa, E., et al. 1991, A&A, 248, 

555 

Young, K. H., Hunter, T. R., Wilner, D. J., et al. 2004, ApJ, 
616, L51 



Menten et al.: Radio continuum observations of IRC+10216 



7 



Zubko, V. G., Mennella, V., Colangeli, L., & Bussoletti, E. 
1996, MNRAS, 282, 1321 



in 
+ 
m 

CD 

o 



CN 
O 

+ 

o 
cr 

c7f 



Julian Date (2440000 + ) 



1; 2.5 



S400 
Fl 



5600 



300 



9000 
\ 



-U3 — ^ 

2 Til = =ii\- 

1.5 - 




(A) 



(/) 



in 



40 


-40 

40 


-40 

40 

-40 

10 



_o. 



i K 











1-- 


1 






f-- 



5 - 




I I I I I I I L 

-0.5 0.5 

hfrored Phase (<i>) 



Fig. 1. Radio flux densities of IRC+10216 from 1991 June 

1 to 1993 June 8 (second panel from top) compared to that 
of the cahbrator 0953+254 (top panel). For both sources 
the X-, U-, and K-band integrated flux densities from 
Gaussian fits are plotted bottom to top. Note the diff'crcnt 
flux density scales. Error bars represent the quadratic sum 
of formal fltting errors and estimates of the errors due to 
absolute flux density calibration (10, 15, and 20% for the 
X-, U-, and K-band data points). In the IRC+10216 panel, 
the solid lines represent the results of the model calcula- 
tions described in the text. The third to fifth panels from 
top show K-, U- and X-band residuals in percent-for each 
epoch the flux density predicted by the model was sub- 
tracted from the measured value and divided by the model 
value. X-, U-, and K-band data points are represented as 
crosses, squares, and circles, respectively. In all panels the 
same time range is shown. The abscissa gives Julian day 

mimbfir in t.ViA inn nnnpl pnrl infrnrfirl iT/-)n^r\A\ Qt.Allar 
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10 100 1000 
Frequency (GHz) 

Fig. 2. Interferometric flux densities of IRC+10216 from 
cm to sub-mm wavelengths. The X-, U-, and K-band cm- 
wave data points are the weighted averages of data from 
six epochs (see Table 1). These error bars are quadratic 
sums of formal uncertainties and absolute calibration er- 
rors, which are assumed to be 10, 15, and 20% for the 
X-, U-, and K-band data, respectively. The error bars of 
the mm and submm wavelength data are discussed in the 
text. The straight line has a slope of 2.0, the blackbody 
value. A least squares fit to all the shown data delivers a 
slope of 1.96 ±0.04. 




1 10^ 10^ 10^ lO"* 

Wavelength [/U.m] 



Fig. 3. Spectral energy distribution for IRC-l-10216. Data 
(squares): J- to M- band from Le Bertre (1987); 12 to 
lOO^im from IRAS; 450Aini from Young et al. (2004) for 2" 
beam; 1.24mm fluxes from Lucas & Guelin (1999), lower 
point for l."6 beam, upper point for total envelope; 88% 
of the total 3mm flux (same authors) come from the inner 
2. "4; 1.3 cm point, this paper. Models: The solid curve is 
the model described in the appendix for maximum lumi- 
nosity, the dotted lines are for minimum light, when the 
stellar luminosity is 2.5 times weaker and the surface tem- 
perature is 1700 K. The dash-dots, which overshoot the 
near IR data points, are flux densities from the star alone. 
The dashed line gives only the stellar flux, after extinction 
by the envelope of Ay — 48 mag. 
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Fig. 4. Normalized model flux density as a function of 
aperture at five wavelengtlis. Tlie star has an assumed 
angular diameter of 0." 075; at a comparable resolution es- 
sentially all the emission comes from the star. Substantial 
dust emission is indicated for beam sizes > 1" at IR wave- 
lengths. 



